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Background:  Accurate  prehospital 
triage  of  trauma  patients  is  difficult,  espe¬ 
cially  in  mass  casualty  situations.  Accord¬ 
ingly,  the  U.S.  Military  has  initiated  a  pro¬ 
gram  directed  toward  improving 
noninvasive  prehospital  triage  algorithms 
based  on  available  physiologic  data.  The 
purpose  of  this  study  was  to  assess  heart 
rate  variability  and  its  association  with 
mortality  in  prehospital  trauma  patients. 

Methods:  Trauma  patients  without 
significant  head  injury  requiring  helicop¬ 
ter  transport  were  identified  from  a  ret¬ 
rospective  research  database.  An  equal 
number,  unmatched  sample  of  patients 
who  lived  were  compared  with  those  who 
died  (n  =  15  per  group).  All  patients  were 
transported  to  a  single  Level  I  urban 
trauma  center.  The  primary  independent 
variable  was  mortality.  Patients  with  Ab¬ 
breviated  Injury  Scale  head  scores  >2 
were  excluded  from  the  analysis,  so  that 
the  effects  seen  were  based  on  hemor¬ 
rhagic  shock.  Age,  sex,  Glasgow  Coma 


Scale  score  (GCS),  blood  pressure,  pulse 
pressure,  pulse,  intubation  rate,  Sp02 
mechanism  of  injury,  transport  time,  and 
time  of  death  after  admission  were  re¬ 
corded.  R-waves  from  the  first  available 
120  seconds  of  usable  data  were  detected 
from  normal  electrocardiograms  and 
heart  rate  variability  was  assessed. 

Results:  Patients  who  died  demon¬ 
strated  a  lower  GCS  (7.9  ±  1.4  versus 
14.4  ±  0.2;  p  =  0.0001)  and  higher  intu¬ 
bation  rate  (53%  of  patients  who  died  ver¬ 
sus  0%  patients  who  lived).  Pulse  rate, 
arterial  pressure,  and  Sp02  were  not  dis¬ 
tinguishable  statistically  between  groups 
(p  =  0.08),  but  pulse  pressure  was  lower  in 
patients  who  died  (39  ±  3  versus  50  ±  2 
mm  Hg;  p  =  0.01).  Compared  with  pa¬ 
tients  who  lived,  those  who  died  had  lower 
normalized  low-frequency  (LF)  power 
(42  ±  6  versus  62  ±  4  LFnu;  p  =  0.009), 
higher  high-frequency  (HF)  power  (42  ±  3 
versus  32  ±  3  HFnu;  p  =  0.04)  and  higher 
HF-to-LF  ratio  (144  ±  30  versus  62  ± 


llnu;  p  =  0.01).  With  absolute  HF/LF  ad¬ 
justed  for  GCS,  the  intergroup  variance 
accounted  for  by  HF/LF  was  reduced  to 

6%  (p  =  0.16). 

Conclusions:  Analysis  of  heart  rate 
variability  provides  insight  into  adequacy 
of  autonomic  compensation  to  severe 
trauma.  In  our  cohort  of  trauma  patients, 
low  pulse  pressures  coupled  with  rela¬ 
tively  higher  parasympathetic  than  sym¬ 
pathetic  modulation  characterized  and 
separated  patients  who  died  versus  pa¬ 
tients  who  survived  traumatic  injuries 
when  standard  physiologic  measurements 
are  not  different.  These  data  do  not  sug¬ 
gest  advantages  of  heart  rate  variability 
analysis  over  GCS  scores,  but  suggest  fu¬ 
ture  possibilities  for  remote  noninvasive 
triage  of  casualties  when  GCS  scores  are 
unattainable. 
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Traumatic  injury  is  the  leading  cause  of  death  from  ages  1 
to  44  years,  and  approximately  40%  of  patients  suffering 
fatal  traumatic  injuries  die  before  they  reach  a  hospital.1 
In  the  military  environment,  this  number  can  be  up  to  90%. 2 
In  civilian  trauma  patients,  arterial  pressure,  heart  rate,  arte¬ 
rial  hemoglobin  saturation,  and  mentation  are  routinely  mon- 
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itored  during  transport;  when  abnormal,  they  clearly  are  re¬ 
lated  to  mortality  and  provoke  rapid  evacuation  and 
immediate  interventions.  In  the  military  or  civilian  mass 
casualty  event,  this  level  of  monitoring  may  be  difficult. 
Additionally,  these  physiologic  changes  are  relatively  late 
secondary  (decompensated  physiologic  responses),  rather 
than  primary,  manifestations  of  hemorrhagic  shock,  and  as 
such  may  not  provide  the  first  responder  with  adequate  in¬ 
formation  regarding  triage  categories,  evacuation  priority, 
and  required  interventions. 

Autonomic  neural  modulation  of  the  circulation  is  the 
primary  compensatory  response  to  traumatic  hemorrhage, 
and  patients  with  adequate  autonomic  responses  are  better 
able  to  maintain  peripheral  resistance  and  arterial  pressure 
compared  with  those  patients  whose  autonomic  responses  are 
not  adequate.  We  have  shown  that  estimates  of  autonomic 
activity  change  as  direct  inverse  functions  of  the  magnitude 
of  negative  pressure  applied  to  the  lower  body  of  human 
volunteers.3  Lower-body  negative  pressure  is  useful  as  an 
experimental  procedure  to  induce  controlled,  progressive  re¬ 
ductions  of  central  blood  volume  and  resulting  cardiovascular 
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compensations  similar  to  actual  hemorrhage.4  Thus,  it  is 
reasonable  to  propose  that  monitoring  of  autonomic  function 
in  trauma  patients  may  provide  early  information  to  the  first 
responder  not  contained  in  other  measurements  of  hemody¬ 
namic  stability. 

Power  spectral  analysis  has  emerged  as  a  useful  nonin- 
vasive  tool  using  simple  electrocardiogram  (ECG)  recordings 
for  quantifying  relative  parasympathetic  and  sympathetic  ac¬ 
tivities  by  separating  R-R  interval  oscillations  into  their  high- 
frequency  (HF;  0.15  to  0.4  Hz)  and  low-frequency  (LF;  0.05 
to  0.15  Hz)  components.  Analysis  of  heart  rate  variability 
provides  quantitative  estimates  of  autonomic  neural  activity5 
and  may  have  prognostic  value  for  hospitalized  trauma 
patients.6-8 

For  example,  in  a  large  cohort  of  both  cardiac  and  non¬ 
cardiac  patients  (n  =  740),  Winchell  and  Hoyt7  estimated 
relative  parasympathetic  and  sympathetic  activities  by  apply¬ 
ing  power  spectral  analysis  to  cardiac  interbeat  intervals. 
They  concluded  that  survival  was  associated  with  low  para¬ 
sympathetic  and  high  sympathetic  predominance,  and  death 
was  associated  with  high  parasympathetic  and  low  sympa¬ 
thetic  predominance.  Heart  rate  variability  analysis  may  char¬ 
acterize  adequacy  of  autonomic  compensation  to  an  acute 
traumatic  event,  but  to  our  knowledge,  heart  rate  variability 
analysis  has  not  been  applied  to  trauma  patients  in  a  prehos¬ 
pital  setting.  The  purpose  of  this  study  was  to  test  the  hy¬ 
pothesis  that  mortality  in  prehospital  trauma  patients  is  asso¬ 
ciated  with  higher  parasympathetic  and  lower  sympathetic 
autonomic  activity  as  estimated  from  analysis  of  heart  rate 
variability.  We  introduce  the  concept  that  estimates  of  auto¬ 
nomic  function  may  allow  for  the  separation  and  classifica¬ 
tion  of  trauma  patients  who  may  appear  normal  utilizing 
standard  physiologic  assessments,9  but  who  will  ultimately 
die  from  their  injuries. 

PATIENTS  AND  METHODS 

This  study  was  approved  by  the  Institutional  Review 
Board  of  the  University  of  Texas  Health  Science  Center, 
Houston,  Texas,  and  the  Research  Review  Committee  of  the 
U.S.  Army  Institute  of  Surgical  Research  and  the  Brooke 
Army  Medical  Center  Institutional  Review  Board,  Fort  Sam 
Houston,  Texas.  Patients  were  enrolled  for  this  study  using 
the  Trauma  Vitals  database  system  developed  by  the  U.S. 
Army  Institute  of  Surgical  Research.  The  system  provides  a 
data  warehousing  capability  for  storing  and  correlating  pre¬ 
hospital  patient  data  from  incident  pickup  until  delivery  via 
Fife  Flight  helicopter  to  Memorial  Hermann  Hospital,  a  re¬ 
gional  Fevel  I  trauma  center  in  Houston.  Fife  Flight  utilizes 
standard  treatment  protocols  to  guide  onboard  interventions. 
From  a  total  patient  population  of  914,  we  identified  93 
patients  who  died  from  their  injuries  ( — 10%).  Of  these  pa¬ 
tients  who  died,  we  excluded  66  patients  because  they  had 
either  incomplete  data,  no  ECG  recorded,  or  had  an  Abbre¬ 
viated  Injury  Score  (AIS)  head  score  >2.  Of  the  remaining  26 
patients  who  met  our  criteria  for  inclusion,  1 1  were  excluded 


Original  Raw  Data 


Fig.  i.  Raw  ECG  data  obtained  directly  from  the  ProPaq  206EL  is 
shown  in  the  upper  panel  (original  raw  data).  The  same  data  are 
shown  in  the  lower  panel  (processed  data)  after  extensive  filtering 
established  stable  baselines  and  removed  noise  spikes  and  move¬ 
ment  artifacts  from  the  signal. 

based  on  the  very  poor  quality  of  their  ECG  signal.  We  were 
therefore  able  to  obtain  analyzable  data  from  1 5  patients  who 
died  from  their  injuries.  We  chose  for  analysis  an  equal 
number  of  15  patients  with  similar  exclusion  criteria  who 
lived.  For  this  exploratory  study,  we  did  not  attempt  to  match 
patient  groups  for  age,  sex,  or  some  other  parameter. 

Data  Collection 

Data  were  collected  using  both  an  automated  data  col¬ 
lection  device  for  capturing  values  generated  by  the  on-board 
vital  signs  monitor  and  through  the  standard  run  sheet  for 
every  patient.  Vital  signs,  Glasgow  Coma  Scale  score  (GCS), 
age,  sex,  mechanism  of  injury,  and  incident  demographics 
were  recorded.  Arterial  pressures  were  measured  automati¬ 
cally  and  recorded  every  3  minutes.  Arterial  pressure  and 
Sp02  (pulse  oximetry)  data  presented  in  this  report  were 
taken  within  3  minutes  of  the  ECG  dataset  used  for  heart  rate 
variability  analysis.  ECG  data  were  collected  with  an  iPAQ 
(Talla-tech  RPDA,  Tallahassee,  FF)  personal  digital  assistant 
interfaced  with  a  ProPaq  206EF  vital  signs  monitor  (Wel- 
chAllyn,  Beaverton,  OR)  and  recorded  with  a  sampling  fre¬ 
quency  of  182  Hz.  All  descriptive  data  including  clinical 
outcomes  were  recorded  prospectively  by  a  single  research 
nurse. 

Heart  Rate  Variability  Analyses 

ECG  data  were  extracted  from  the  Trauma  Vitals  data¬ 
base,  and  then  manually  filtered  and  cleaned  with  customized 
software  to  remove  noise  artifacts.  An  example  of  the  mag¬ 
nitude  of  data  filtering  required  for  one  representative  patient 
is  shown  in  Figure  1 . 

Filtered  and  cleaned  ECG  signals  were  analyzed  with  a 
commercially-available  software  analysis  package 
(WinCPRS,  Absolute  Aliens,  Turku,  Finland).  We  identified 
the  first  2  minutes  of  continuous  data  for  each  patient,  and 
these  sections  were  considered  valid  if  data  were  clean 
enough  to  identify  with  certainty  individual  R  waves  for  each 
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Heart  Rate  Variability  and  Trauma 


Table  1  Patient  Characteristics 


Variable 

Died 

Lived 

F  Ratio 

p  Value 

n 

15 

15 

Age  (years) 

43  ±  2 

35  ±  3 

2.2 

0.15 

Sex 

Male 

13 

13 

Female 

2 

2 

Blunt  injury  (%) 

73 

60 

0.6* 

0.45 

Glasgow  Coma  Scale  score 

7.9  ±  1.4 

14.4  ±  0.2 

19.8 

0.0001 

Time  to  death  (hours) 

153  ±  59 

Transport  time  (minutes) 

15.3  ±  1.8 

12.4  ±  1.5 

1.8 

0.19 

Estimated  time  after  injury  to  data  collection  (minutes) 

42.2  ±  3.6 

49.2  ±  4.1f 

1.5 

0.22 

Intubated  (%) 

53 

0 

10.9* 

0.001 

Values  are  means  ±  SE. 
'Estimated  by  x2  test. 
fn  =  14. 


cardiac  cycle.  Data  sets  including  noisy  ECGs,  ectopic  beats, 
or  “spikes”  that  could  not  be  identified  clearly  by  their  posi¬ 
tion  in  the  cardiac  cycle  were  excluded  from  analysis  and 
these  patients  were  not  included  in  the  final  dataset:  a  total  of 
1 1  of  26  patients  who  died  and  met  our  inclusion  criteria  were 
excluded  from  analysis  due  to  these  criteria  for  clean  ECG 
signals. 

Heart  rate  variability  was  assessed  in  the  time  domain  by 
calculating  the  standard  deviation  (SD)  of  average  R-R  in¬ 
tervals  (RRISD),  the  square  root  of  the  mean  squared  differ¬ 
ences  of  successive  R-R  intervals  (RMSSD),  and  the  per¬ 
centage  of  adjacent  normal  R-R  intervals  that  varied  by  at 
least  50  ms  or  more  (pNN50).  Heart  rate  variability  was 
assessed  in  the  frequency  domain  from  R-R  interval  spectral 
power.  To  calculate  power  spectra,  consecutive  R-R  intervals 
were  made  equidistant  by  spline  interpolating  and  resampling 
at  5  Hz.  Data  then  were  passed  through  a  low-pass  impulse 
response  filter  with  a  cutoff  frequency  of  0.5  Hz.  Data  sets 
were  analyzed  in  the  frequency  domain  using  a  Fourier  trans¬ 
form  with  a  Hanning  window  as  described  previously.9-13 
Heart  rate  variability  was  quantified  by  calculating  the  total 
integrated  area  under  the  power  spectrum  (Tot:  0.05-0.4  Hz). 
Integrated  areas  were  also  separated  into  HF  (0.15-0.4  Hz) 
and  LF  (0.05-0.15  Hz)  bands.  To  compare  power  spectra 
more  accurately  from  different  patients  who  may  vary  widely 
in  total  power,  we  normalized  our  data  in  the  FF  (LFnu)  and 
HF  (HFnu)  ranges  by  dividing  integrated  LF  and  HF  spectra 
by  the  total  power  (minus  oscillations  occurring  below  0.05 
Hz)  and  multiplying  this  value  by  100. 1415  We  took  the 
quotient  of  HF/LF  as  a  global  index  of  changes  in  autonomic 
balance.15 

Statistical  Analysis 

We  analyzed  our  data  with  commercial  statistical  soft¬ 
ware  (SAS  Institute,  Cary,  NC).  The  primary  independent 
variable  for  univariate  analyses  was  mortality.  Differences 
between  the  means  of  each  dependent  variable  were  tested 
with  a  one-way  ANOVA  with  comparisons  between  groups 


(patients  who  died  versus  patients  who  lived).  We  also  ana¬ 
lyzed  our  data  further  using  logistic  regression.  For  the  mul¬ 
tivariate  phase,  the  lived/died  grouping  was  modeled  as  the 
dependent  effect.  The  variable  HF/LF  was  found  to  contain 
all  of  the  relevant  information  from  each  of  the  redundant 
frequency-domain  calculations  performed  for  the  univariate 
analyses.  Logistic  regression  was  used  to  test  for  the  unique 
effect  of  HF/LF  after  controlling  for  each  of  three  covariates 
(age,  sex,  and  GCS  score).  To  aid  in  the  interpretation  of  the 
results,  adjusted  group  means  were  calculated  via  inverse 
prediction  of  the  logistic  model.  Exact  p  values  were  calcu¬ 
lated  for  all  statistical  tests  and  reflect  the  probability  of 
obtaining  the  observed  or  greater  effect  given  only  random 
departure  from  the  assumption  of  no  effects.  All  data  are 
expressed  as  means  ±  standard  errors  unless  specified  oth¬ 
erwise. 

RESULTS 

Subject  characteristics,  classifications,  and  descriptive 
data  are  displayed  in  Table  1.  Injury  mechanisms  (blunt  or 
penetrating)  were  similar  between  groups.  Transport  time  and 
the  estimated  time  after  incident  occurrence  until  admission 
to  the  emergency  department  were  not  distinguishable  statis¬ 
tically  between  patient  populations  (Table  1).  Median  time  to 
death  after  injury  for  patients  who  died  was  9.5  hours.  Pa¬ 
tients  who  died  had  lower  GCS  scores  compared  with  patients 
who  lived.  Intubation  rate  was  53%  for  patients  who  died;  no 
patients  who  lived  were  intubated.  No  patient  who  died  had 
an  AIS  head  score  >1. 

Electrocardiograms  obtained  during  transport  on  the  Life 
Flight  helicopter  were,  as  expected,  consistently  noisy  and 
contained  numerous  movement  artifacts.  To  obtain  a  sample 
size  appropriate  for  rigorous  statistical  analysis,  we  deter¬ 
mined  a  priori,  based  on  the  availability  of  clean  ECG  signals 
from  which  R  waves  were  clearly  distinguishable,  that  anal¬ 
yses  would  be  done  on  2-minute  data  sets.  Two  minutes  is 
sufficient  to  interpret  oscillations  occurring  within  the  fre¬ 
quency  band  of  interest  (0.05-0.4  Hz).14 
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T3bl6  2  Time  Domain  Results  for  Patients  Who  Died 
and  Patients  Who  Survived  Traumatic  Injuries  during 
Transport  on  a  Life  Flight  Helicopter 

Variable 

Died 

Lived 

F  Ratio 

p  value 

n 

15 

15 

RRI  (ms) 

605  ±  37 

712  ±  46 

3.3 

0.08 

RRISD  (ms) 

58  ±  14.6 

47  ±  7 

0.45 

0.51 

RMSSD  (ms) 

53  ±  12 

43  ±  7 

0.47 

0.50 

pNN50 (%) 

13  ±  5 

13  ±  4 

0.0 

0.96 

Heart  rate  (bpm) 

105  ±  6 

90  ±  6 

3.2 

0.08 

Sp02 (%) 

90  ±  6* 

98  ±  0.9f 

2.3 

0.14 

Values  are  means  ±  SE. 
*n  =  9. 
tn  =  12. 


Table  3  Frequency  Domain  Results  for  Patients  Who 
Died  and  Patients  Who  Survived  Traumatic  Injuries 
during  Transport  on  a  Life  Flight  Helicopter 


Variable  Died  Lived  F  Ratio  p  Value 


n 

15 

15 

HF  (ms2) 

635 

+ 

434 

457 

+ 

163 

0.15 

0.70 

LF  (ms2) 

618 

+ 

322 

833 

+ 

272 

0.26 

0.61 

Tot  (ms2) 

1253 

+ 

750 

1290 

-+- 

411 

0.0 

0.96 

LF/HF 

240 

+ 

147 

250 

-+- 

46 

0.0 

0.95 

HF/LF 

174 

+ 

32 

65 

-+- 

12 

10.8 

0.003 

HFnu 

42 

+ 

3 

32 

+ 

3 

4.3 

0.04 

LFnu 

42 

HF 

6 

62 

+ 

4 

7.7 

0.009 

HF/LFnu 

146 

± 

30 

62 

+ 

11 

6.9 

0.01 

LF/HFnu 

129 

+ 

28 

250 

-+- 

45 

5.1 

0.03 

RRI,  R-R  interval;  RRISD,  R-R  interval  standard  deviation; 
RMSSD,  the  square  root  of  the  mean  squared  differences  of  succes¬ 
sive  R-R  intervals;  pNN50,  the  percentage  of  adjacent  normal  R-R 
intervals  that  vary  by  at  least  50  ms  or  more;  Sp02,  arterial  oxygen 
saturation. 


Numerical  results  derived  from  time-domain  analyses 
are  shown  in  Table  2.  Heart  rates  and  R-R  intervals,  as  well 
as  global  measures  of  heart  rate  variability  including  RRISD, 
RMSSD,  and  pNN50  were  indistinguishable  statistically  be¬ 
tween  the  two  patient  groups.  Sp02  values  were  unstable  and 
inconsistent  for  nine  patients.  Values  were  determined  to  be 
stable  and  reliable  when  pulse  rates  from  the  pulse  oximeter 
matched  pulse  rates  from  the  radial  pulse.  These  Sp02  values 
were  taken  within  3  minutes  of  the  ECG  data  used  to  assess 
heart  rate  variability.  With  these  criteria,  we  analyzed  Sp02 
from  nine  patients  who  died  and  12  patients  who  lived,  and 
found  that  Sp02  was  similar  between  groups.  Patient  groups 
had  similar  systolic,  diastolic,  and  mean  arterial  pressures, 
but  pulse  pressures  were  lower  (39  ±  3  versus  50  ±  2  mm 
Hg;  p  =  0.01)  in  patients  who  died  compared  with  patients 
who  lived  (Fig.  2). 

Table  3  shows  results  from  frequency-domain  analyses. 
Total  power,  HF  and  LF  power  were  not  different  statistically 
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Fig.  2.  Arterial  pressures  and  pulse  pressures  are  shown  for  13 


Values  are  means  ±  SE.  HF,  R-R  interval  spectral  power  at  the 
high  frequency  (0.15-0.4  Hz);  LF,  R-R  interval  spectral  power  at  the 
low  frequency  (0.05-0.15  Hz);  Tot,  total  R-R  interval  spectral  power 
(0.05-0.4  Hz);  LF/HF,  ratio  indicating  relative  sympathetic  predomi¬ 
nance;  HF/LF,  ratio  indicating  relative  parasympathetic  predomi¬ 
nance;  HFnu,  normalized  high  frequency  R-R  interval  spectral  power; 
LFnu,  normalized  low  frequency  R-R  interval  spectral  power;  HF/LFnu, 
normalized  ratio  indicating  relative  parasympathetic  predominance; 
LF/HFnu,  normalized  ratio  indicating  relative  sympathetic  predomi¬ 
nance. 

between  the  two  groups  when  data  were  expressed  as  abso¬ 
lute  values.  However,  it  is  recommended  that  comparisons 
between  spectral  components  also  be  made  on  data  that  have 
been  normalized  for  the  absolute  value  of  the  total  variance  of 
the  signal  (analogous  to  total  power).15  HFnu  and  the  HF/LFnu 
ratio  (including  the  ratio  expressed  in  absolute  units)  was 
higher,  LFnu  and  the  LF/HF11U  ratio  was  lower  in  patients  who 
died  from  their  injuries  compared  with  patients  who  lived 
(Table  3). 

Figure  3  shows  average  spectral  power  for  all  30  patients 


N 

I 

N* 

(f) 

E 

Q 

CO 

CL 

rE 

QL 


0.05  0.10  0.15  0.20  0.25  0.30  0.35  0.40 

frequency,  Hz 

Fig.  3.  Average  R-R  inten’al  power  spectral  density  (RRI  PSD)  is 
shown  for  15  patients  who  lived  (solid  line )  and  15  patients  who  died 
( broken  line  with  red  shading)  traumatic  injuries  requiring  trans¬ 
port  to  a  level  one  trauma  center;  vertical  dotted  line  denotes  the 


patients  who  died  and  15  patients  who  survived  traumatic  injuries 


demarcation  of  low-frequency  (0.05-0.15  Hz)  and  high-frequency 


requiring  transport  to  a  level  one  trauma  center.  *p  =  0.001.  (0.15-0.4  Hz)  spectral  bands. 
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Table  4  Multivariate  Summary 


Adjusted 

Effect  Model  ^rolJP  Likelihood  Ratio 

Tested  Covariate(s)  Means  (x^hf/lf ) 


Died  Lived 


HF/LF 

— 

174 

65 

9.96 

0.0016 

HF/LF 

Age 

171 

68 

5.19 

0.0227 

HF/LF 

RRI 

167 

72 

7.06 

0.007 

HF/LF 

GCS 

152 

87 

1.96 

0.1619 

HF/LF 

Age,  RRI,  GCS 

132 

108 

0.43 

0.487 

Multivariate  model  with  between  group  difference  in  HF/LF  ad¬ 
justed  for  variables  that  also  indicated  between  group  differences. 

HF,  R-R  interval  spectral  power  at  the  high  frequency  (0.15-0.4 
Hz);  LF,  R-R  interval  spectral  power  at  the  low  frequency  (0.05-0.15 
Hz);  RRI,  interbeat  R-R  interval;  GCS,  Glasgow  Coma  Scale  score. 

separated  into  those  who  lived  and  those  who  died.  Figure  3 
shows  graphically  the  differences  in  R-R  interval  oscillations 
at  low  frequencies  around  0.1  Hz  and  at  high  frequencies 
around  0.3  Hz. 

The  univariate  analysis  only  suggests  possible  differ¬ 
ences  that  may  exist  between  subjects  who  die  and  those  who 
survive.  Such  analysis  does  not  account  for  the  redundancies 
that  exist  between  variables  and  the  confounds  that  exist 
between  groups.  We  determined  that  from  all  of  the  frequen¬ 
cy-domain  variables,  HF/LF  in  absolute  units  accounts  for  the 
greatest  amount  of  between  group  variance.  Table  4  presents 
the  results  of  a  multivariate  model  in  which  the  between 
group  difference  seen  in  HF/LF  was  adjusted  for  variables 
that  also  indicated  between  group  differences  (Tables  1  and 
2).  These  covariates  included  age,  GCS,  and  RRI.  When 
HF/LF  is  adjusted  for  age  or  RRI,  there  is  a  slight  reduction 
in  HF/LF  group  means  compared  with  the  HF/LF  univariate 
results.  However,  the  HF/LF  effect  still  explains  approxi¬ 
mately  20%  of  the  between  group  variance  when  adjusted  for 
either  age  or  RRI.  When  HF/LF  is  adjusted  for  GCS,  the 
between  group  variance  accounted  for  by  HF/LF  is  reduced  to 
6%  and  the  p  value  associated  with  the  HF/LF  effect  drops  to 
0.1619.  When  all  three  covariates  are  entered  into  the  model 
simultaneously,  the  variance  accounted  for  by  the  HF/LF 
effect  and  the  associated  p  values  are  1%  and  0.4877,  respec¬ 
tively.  These  results  are  shown  in  Table  4. 

DISCUSSION 

We  assessed  heart  rate  variability  from  electrocardio¬ 
grams  recorded  in  the  prehospital  environment  from  trauma 
patients  who  died  to  assess  whether  mortality  is  associated 
with  estimates  of  autonomic  function.  The  primary  new  find¬ 
ings  are:  1)  analysis  of  heart  rate  variability  in  prehospital 
trauma  patients  provides  insight  into  adequacy  of  autonomic 
compensation  and  survival  after  trauma;  and  2)  patients  who 
survive  their  injuries  have  higher  pulse  pressures  and  main¬ 
tain  high  sympathetic  modulation,  whereas  patients  who  die 
have  low  pulse  pressures,  withdraw  their  sympathetic  tone 
and  display  inappropriate  parasympathetic  neural  predomi¬ 


nance  as  estimated  from  analysis  of  heart  rate  variability. 
However,  heart  rate  variability  analysis  contributes  no  addi¬ 
tional  predictive  value  in  a  multivariate  model  when  it  is 
covaried  with  GCS  scores.  Our  results  suggest  that  autonomic 
balance  and  pulse  pressure  are  associated  with  mortality,  and 
may  provide  important  diagnostic  and  prognostic  information 
to  assist  the  first  responder  in  on-site  triage,  remote  triage, 
and  evacuation  priority  in  the  absence  of  obtainable  GCS 
scores. 

Oscillations  in  heart-period  intervals  in  humans  separate 
into  HF  (0.15-0.4  Hz)  and  LF  (0.05-0.15  Hz)  bands.14  R-R 
interval  oscillations  at  HF  are  abolished  with  vagal 
blockade,16  and  so  changes  in  R-R  interval  oscillations  at  HF 
represent  primarily  parasympathetic  activity.  Oscillations  oc¬ 
curring  at  frequencies  lower  than  respiration  are  blunted  with 
vagal  blockade  but  abolished  with  the  combination  of  vagal 
and  sympathetic  blockade,16  and  so  LF  R-R  interval  oscilla¬ 
tions  are  modulated  importantly  by  both  parasympathetic  and 
sympathetic  activities.  Because  LF  R-R  interval  rhythms  are 
modulated  by  both  branches  of  the  autonomic  nervous  sys¬ 
tem,  a  vigorous  debate  has  developed  over  the  concept  of 
sympatho vagal  balance  (the  ratio  of  LF/HF).  Eckberg17  has 
argued  that  “calculations  of  sympathovagal  balance  may  ob¬ 
scure  rather  than  illuminate  human  physiology  and  patho¬ 
physiology,”  and  Malliani15  maintains  that  an  indirect  quan¬ 
tification  of  the  relative  contributions  of  sympathetic  and 
parasympathetic  modulation  of  sinus  node  rate  may  be  ob¬ 
tained  by  calculating  the  ratios  of  LF  and  HF  R-R  interval 
spectral  powers.  Strong  cases  have  been  made  both  for15  and 
against17  the  concept  of  sympathovagal  balance,  but  when 
both  LF  and  HF  R-R  interval  powers  are  normalized  to 
control  for  changes  in  total  signal  variance,  LFnu  increases 
consistently,  and  HFnu  decreases  consistently  with  maneuvers 
that  increase  sympathetic  and  decrease  parasympathetic  neu¬ 
ral  activities.18  In  trauma  care,  the  argument  over  whether 
spectral  components  of  cardiac  interbeat  intervals  track  ab¬ 
solute  activities  of  autonomic  nerves  is  moot  if  interpretations 
of  such  spectral  components  give  insights  into  patient  out¬ 
come. 

Compensations  for  traumatic  injuries  with  consequent 
reductions  of  central  blood  volume  include  increased  heart 
rate  and  peripheral  vascular  resistance.  These  autonomic 
compensations  are  mediated  by  withdrawal  of  parasympa¬ 
thetic  and  activation  of  sympathetic  efferent  neural  traffic  to 
the  heart  and  vasculature.  Controlled  experimental  hemor¬ 
rhage  in  animals  has  revealed  two  primary  phases  of  auto¬ 
nomic  compensation:  1)  sympathoexcitation,  resulting  in 
maintenance  of  arterial  pressure  and  central  organ  perfusion 
by  vasoconstriction;  and  2)  pronounced  sympathoinhibition, 
resulting  in  vasodilation,  bradycardia,  and  eventual  cardio¬ 
vascular  collapse.19  Such  phases  are  consistent  with  the  ob¬ 
servations  of  Barcroft  and  coworkers,20  who  demonstrated  in 
a  seminal  study  that  circulatory  collapse  due  to  hemorrhage 
occurs  in  conjunction  with  sudden  drops  of  peripheral  resis- 
tance.  We"  and  others"  ’  have  reported  previously  that 
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during  induced  hypotension,  peripheral  resistance  falls  con¬ 
sequent  to  abrupt  withdraw  of  sympathetic  traffic.  In  humans 
undergoing  simulated  hemorrhage  with  lower  body  negative 
pressure21'23  or  induced  hypovolemia  with  passive  head-up 
tilt  and  nitroprusside  infusion,22  arterial  pressures  were  main¬ 
tained  through  sympathoexcitation  up  to  the  point  where 
further  increases  in  sympathetic  activity  were  not  possible. 
When  subjects  reached  their  limits  for  sympathetic  activation, 
they  experienced  profound  acute  reductions  of  arterial  pres¬ 
sure  in  conjunction  with  total  withdrawal  of  sympathetic 
traffic  (measured  directly  from  the  peroneal  nerve  with 
microneurography).2 1-2  Results  from  such  controlled  labo¬ 
ratory  experiments  suggest  that  maintenance  of  robust  sym¬ 
pathetic  nerve  traffic  is  a  prerequisite  for  maintaining  com¬ 
pensated  hemodynamics  during  hemorrhage. 

However,  it  is  not  feasible  to  record  directly  and  contin¬ 
uously  from  a  sympathetic  nerve  in  patients  during  either 
transport  to  a  hospital  or  during  prolonged  observation  in  a 
hospital.  Moreover,  direct  measurements  from  the  vagus 
nerve  have  not  been  made  in  humans,  but  it  may  be  possible 
to  estimate  relative  parasympathetic  and  sympathetic  pre¬ 
dominance  non-invasively  by  analyzing  heart  rate  variability. 
Low  total  power  (an  estimate  of  overall  autonomic  tone),  and 
high  HF/LF  (an  estimate  of  parasympathetic  predominance), 
were  found  to  be  associated  with  mortality  in  742  patients 
monitored  over  the  course  of  6  months  in  a  surgical  ICU.7 
Based  on  these  results  obtained  primarily  from  stable  pa¬ 
tients,  we  hypothesized  that  power  spectral  analysis  might 
also  be  of  benefit  in  evaluating  patients  and  providing  prog¬ 
noses  in  a  prehospital  scenario  where  patients  may  not  be 
stable.  Data  presented  in  Table  2  highlight  the  lack  of  differ¬ 
ence  in  heart  rates,  Sp02  (in  a  small  number  of  patients)  or 
global  measures  of  heart  rate  variability  calculated  in  the  time 
domain.  The  primary  benefit  of  frequency  domain  over  time 
domain  analysis  is  the  capacity  to  separate  variabilities  into 
frequency  components  that  are  related  to  physiologic  pro¬ 
cesses.  Normalized  frequency  bands  and  the  autonomic  ratios 
suggest  significantly  greater  parasympathetic  and  lesser  sym¬ 
pathetic  predominance  in  patients  who  died  versus  patients 
who  survived  (Table  3).  Our  results  are  consistent  with  pre¬ 
vious  reports  of  sympathetic  withdrawal  and  associated  hy¬ 
potension  during  simulated  hemorrhage2123  and  induced 
hypotension,22  and  are  in  agreement  with  the  suggestion  that 
mortal  outcomes  in  ICU  patients  are  associated  with  high 
parasympathetic  and  low  sympathetic  neural  tone.6'7 

Patients  who  died  were  clearly  identified  by  the  flight 
medics  as  critically  ill,  as  evidenced  by  an  increased  intuba¬ 
tion  rate  and  decreased  GCS.  These  experienced  providers 
recognized  that  these  casualties  were  severely  injured  and 
required  rapid  evacuation,  coupled  with  immediate  live  sav¬ 
ing  interventions.  It  is  commonly  recognized  that  young 
trauma  patients  can  compensate  for  an  extended  period  of 
time  before  rapidly  decompensating.  The  inexperienced  pro¬ 
vider,  or  one  caring  for  overwhelming  mass  casualty  patients, 
may  not  recognize  the  precarious  patient  with  a  narrowed 


pulse  pressure  that  is  moments  away  from  cardiovascular 
collapse.  The  value  of  the  results  presented  in  this  report  is  in 
the  potential  to  deliver  early  indicators  of  injury  severity, 
perhaps  even  from  a  remote  location  before  the  arrival  of  a 
first  responder.  Heart  rates,  Sp02,  and  blood  pressures  are 
commonly  recorded  by  the  medical  monitors  that  are  placed 
in  emergency  vehicles  and  aircraft.  However,  our  data  sug¬ 
gest  that  these  vital  signs  fail  to  provide  a  first  responder  with 
early  adequate  information  as  to  whether  a  bleeding  patient  is 
progressing  toward  cardiovascular  collapse  that  could  lead, 
ultimately,  to  hemorrhagic  shock  and  death.24 

In  the  present  study,  heart  rates,  Sp02,  and  arterial  pres¬ 
sures  were  not  different  between  the  two  patient  populations 
but  pulse  pressures  were  lower  in  patients  who  died  compared 
with  patients  who  lived.  We  speculate,  but  cannot  prove,  that 
changes  of  pulse  pressure  track  changes  of  central  blood 
volume.  In  a  recent  study,  Leonetti  et  al.25  found  that  both 
stroke  volume  and  pulse  pressure  correlated  closely  with 
blood  loss  after  phlebotomy  therapy  of  approximately  350 
mL.  Although  the  speculation  that  changes  of  pulse  pressure 
track  changes  of  blood  volume  requires  confirmation  with 
animal  models  of  controlled  hemorrhage,  we  suggest  that 
automated  calculation  of  pulse  pressure  rather  than  presenting 
a  mean  arterial  pressure  may  provide  the  first  responder  with 
information  on  primary  volume  status  after  injury  rather  than 
information  pertaining  to  a  secondary  response  to  hemor¬ 
rhage. 

Data  presented  in  this  report,  and  those  of  others,6-8  also 
suggest  that  important  information  regarding  the  integrity  and 
adequacy  of  autonomic  neural  function  may  be  contained  in 
components  of  heart  rate  variability.  However,  when  esti¬ 
mates  of  autonomic  function  are  adjusted  for  injury  severity 
as  quantified  by  GCS  scores,  such  estimates  may  not  provide 
the  first  responder  with  any  additional  useful  information 
(Table  4).  We  propose  that  the  utility  of  using  frequency- 
domain  analysis  in  prehospital  trauma  care  may  be  specific 
for  remote  triage  or  other  conditions  where  visual  inspection 
or  palpation  of  a  casualty  may  not  be  possible.  We  also  stress 
that  there  is  at  least  one  limitation  to  using  heart  rate  vari¬ 
ability  analysis  as  a  diagnostic  tool  in  the  prehospital  patient, 
and  that  is  the  inability  to  obtain  real-time  spectral  estimates 
from  clean  ECG  signals.  In  their  studies  of  patients  with 
severe  brain  injuries,  Baillard  et  al.8  collected  ECGs  under 
stable  conditions  in  an  ICU,  insuring  clean  signals  with 
readily  identifiable  R-wave  fiducial  points.  It  was  possible 
then  to  develop  algorithms  to  provide  R-R  interval  spectral 
components  in  real  time,  a  prerequisite  for  monitoring  and 
interpreting  changes  in  heart  rate  variability  from  some  stable 
baseline  value.8  In  the  present  study,  ECGs  collected  from  the 
ProPaq  206EL  onboard  Life  Flight  helicopters  were,  as  ex¬ 
pected,  consistently  noisy  and  required  extensive  manual 
filtering  and  postprocessing  before  we  could  determine  with 
confidence  R-wave  fiducial  points  (Fig.  1).  To  investigate 
further  the  possibility  of  using  frequency  domain  analysis  of 
R-R  intervals  as  a  diagnostic  and  prognostic  tool  for  the  first 
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responder  in  real  time,  portable  field  vital  sign  monitoring 
devices  must  output  cleaner  ECG  signals  with  readily  iden¬ 
tifiable  R-waves  and  have  built  in  data  analysis  capability. 

A  second  limitation  to  the  current  study  is  the  small 
number  of  patients  compared.  We  recognize  and  stress  that 
clinical  decisions  based  on  heart  rate  variability  analysis 
should  derive  from  robust  samples  of  patients  with  appropri¬ 
ate  statistical  controls.  Until  limitations  associated  with  tech¬ 
nical  issues  and  small  sample  sizes  are  addressed,  our  pre¬ 
liminary  but  unique  observations  only  suggest  at  this  time  a 
reasonable  hypothesis  that  should  be  explored  further.  Re¬ 
sults  of  logistic  regression  analysis  showed  that  frequency 
domain  results  are  confounded  by  GCS  scores.  In  future 
experiments,  survivors  should  be  compared  with  nonsurvi¬ 
vors  for  patients  who  have  the  same  GCS  score.  Matching 
patients  for  GCS  scores  is  necessary  to  determine  whether 
heart  rate  variability  analysis  is  predictive  of,  rather  than 
simply  associated  with,  mortality. 

Based  on  frequency  domain  analysis  of  R-R  intervals  from 
data  collected  en  route  to  a  hospital,  trauma  patients  who  die 
from  hemorrhagic  shock  increase  their  parasympathetic  and  de¬ 
crease  their  sympathetic  predominance  compared  with  patients 
who  survive  their  injuries.  Although  we  concede  that  these 
results  must  be  confirmed  with  larger  patient  populations,  our 
results  are  consistent  with  reports  that  hemodynamic  instability 
during  experimentally-induced  hypovolemia  is  associated  with 
sympathetic  neural  withdrawal,  and  also  with  reports  that  mor¬ 
tality  in  ICU  patients  is  associated  with  high  parasympathetic 
and  low  sympathetic  autonomic  tone.  Such  autonomic  responses 
occur  in  conjunction  with  significantly  lower  pulse  pressures  in 
patients  who  die,  suggesting  that  pulse  pressure  may  provide  the 
first  responder  with  a  more  accurate  assessment  of  blood  loss 
than  arterial  pressure  alone.  In  the  absence  of  obtainable  GCS 
scores,  analysis  of  heart  rate  variability  represents  a  potential 
new  approach  to  remote  triage.  Future  devices  incorporating 
real-time  frequency  domain  analysis  of  R-R  intervals  and  con¬ 
tinuous  estimates  of  pulse  pressure  could  be  constructed  that 
would  assist  the  first  responder  in  identifying  those  casualties 
that  require  triage,  evacuation,  and  intervention  priority. 
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EDITORIAL  COMMENT 

In  this  issue  of  the  Journal  of  Trauma,  Cooke  and  col¬ 
leagues  report  their  preliminary  observations  of  prehospital 
physiologic  variability  in  15  seriously  injured  patients  who 
died  and  15  matched  survivors.  They  performed  retrospective 
analysis  of  continuous  electrocardiographic  traces  that  were 
acquired  and  stored  in  transport  monitors.  The  authors  as¬ 
sessed  generally  accepted  time-domain  (such  as  standard  de¬ 
viation)  and  frequency-domain  (such  as  high-  and  low-fre¬ 
quency  power)  statistics.  The  major  finding  of  the  study 
related  to  different  high-  and  low-frequency  power  ratios  in 
patients  who  died  compared  with  patients  who  survived. 

What  is  important  about  this  study  is  its  attempt  to 
describe  the  dynamics  of  what  Samuel  V.  Gross  famously 
described  as  the  “rude  unhinging  of  the  machinery  of  life”  in 
the  prehospital  setting  of  traumatic  shock.  In  contrast  to 
complementary  reports  that  focus  on  the  cardiac  state  after 
substantial  resuscitation  (for  example,  in  the  intensive  care 
unit),  the  present  study  explores  endogenous  responses  dur¬ 
ing  the  earliest  phases  of  emergency  care. 

The  data,  which  reflect  comparatively  small  numbers  of 
patients,  suggest  that  the  ratio  of  high-frequency  (HF)  to 
low-frequency  (LF)  power  is  different  among  patients  who 
die  compared  with  patients  who  survive.  High  HF/LF  ratios 
were  associated  with  mortality.  The  authors  follow  conven¬ 
tional  wisdom,  interpreting  a  high  HF/LF  ratio  as  an  auto¬ 
nomic  state  in  which  parasympathetic  activity  predominates 
over  sympathetic  tone.  In  contrast,  survivors  had  a  lower 
HF/LF  ratio,  which  is  commonly  interpreted  as  a  robust 
sympathetic  response.  Put  differently,  the  patients  who  died 
appeared  to  have  exhausted  the  sympathetic  arm  of  the  au¬ 
tonomic  nervous  system.  Such  sympathetic  exhaustion  has 
been  suggested  by  others  to  represent  a  preterminal  condition. 
Indeed,  the  authors’  own  analysis  suggests  that  a  high  HF/LF 


ratio  adds  little  to  a  low  Glasgow  Coma  Scale  score,  the  latter 
also  predicting  death. 

What,  then,  is  the  potential  clinical  value  of  HF/LF 
measures?  First,  they  can  be  made  continuously.  Second,  they 
do  not  require  a  trained  observer  for  assessment,  nor  should 
they  be  subject  to  interrater  variability.  Third,  with  modern 
telemetry,  they  could  be  obtained  remotely  and  simulta¬ 
neously  in  large  numbers  of  patients  or  potential  patients, 
such  as  warfighters  on  the  battlefield.  In  mass  casualty  situ¬ 
ations,  such  measures  could  conceivably  be  used  to  triage 
victims. 

Implementation  and  deployment  of  technology  for  such 
measures  as  clinical  routine  would  have  to  overcome  signif¬ 
icant  technical  hurdles  (the  power  spectral  analysis  that  is 
used  to  generate  the  HF/LF  measure  is  sensitive  to  artifact) 
and  demonstrate  its  value  in  much  larger  studies.  There  are 
challenges  related  to  the  influence  of  chronic  illnesses  and 
medications  on  the  range  of  normal  values  for  HF/LF  ratios. 
There  are  no  population  based  studies  of  stress-related  vari¬ 
ation  in  HF/LF  that  yet  inform  us  as  to  the  presence  or 
absence  of  heritable  factors  that  might  affect  HF/LF.  Put 
simply,  there  is  a  great  deal  to  learn  about  how  this  metric 
might  perform  as  a  predictor  in  a  larger  and  more  heteroge¬ 
neous  population. 

Perhaps  the  most  intriguing  aspect  of  the  data  is  the  way 
it  informs  the  basic  physiologic  sciences.  The  “rude  unhing¬ 
ing”  now  appears  to  be  nonrandom:  Cooke  and  colleagues’ 
data  hint  that  death  has  a  specific  dynamic,  a  dynamic  with  a 
time  signature  that  can  be  extracted  from  conventional  phys¬ 
iologic  monitoring.  If  this  notion  is  correct — if  the  beginning 
of  the  “death  spiral”  can  be  detected  before  it  is  apparent  to 
an  observer  or  manifests  in  vital  signs — then  decompensation 
should  be  identifiable  (and  hopefully  remediable)  before 
complete  physiologic  collapse  becomes  both  inevitable  and 
irreversible. 
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